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Motivated by the wide ranging experimental results on the conductivity of DNA, we have investi- 
gated extraordinary configurations and chemical environments in which DNA might become a true 
molecular wire, perticularly from enhanced electronic overlaps or from small activation energies. 
In particular, we examine A- vs B-DNA, the ribbon-like structures proposed to arise from molec- 
ular stretching, the potential role of counterions in hole doping the DNA orbitals, the possibility 
of backbone conduction, and the effects of water. We find that small activation gaps observed in 
conductivity experiments may arise in the presence of water and counter ions. We further discuss 
the role of harmonic vibration and twisting motion on electron tight binding matrix elements using 
ab initio density functional theory and model Koster-Slater theory calculations. We find that par- 
tial cancellation between ppa and ppir interaction of pz orbitals on adjacent base pairs, along with 
destructive interference of phase factors are needed to explain the weak conductance of A-DNA. 
Our results lead also to a physical interpretation of the angular dependence of inter-base pair tight 
binding matrix elements. Furthermore, we estimate Franck-Condon factors, reorganization energies 
and nuclear frequencies essential for charge transfer rates, and find our estimated hole transfer rates 
between base pairs to be in excellent agreement with recent picosecond dynamics data. 



I. INTRODUCTION 

In addition to DNA's fundamental role in genetics, it 
is now also a potential candidate for nano-electronic de- 
vices. The highly specific binding between single strands 
of DNA and its self-assembly property open a whole new 
approach to single-molecule electronics. However, its in- 
trinsic conductance properties remain highly controver- 
sial. As early as 1962, Elcy and Spivey suggested that 
TT— TT interactions of stacked base pairs in double stranded 
DNA could lead to conducting behavior[l]. Recently, 
there have been several experimental studies of the DNA 
conductance, leading to a variety of results. These have 
ranged from wide-gap insulating behavior to proximity 
induced superconductivity. 

Although double stranded DNA has similarities to van 
der Waals stacked aromatic conducting crystals {e.g. the 
DNA base pair spacing is similar to the preferred axis 
lattice spacing of aromatic crystals) it also has essential 
differences from these and conventional conductors. Un- 
like crystals, the DNA is not a periodic system. The 
largest ionization potential difference between two iso- 
lated bases is about 0.6 eV between guanine and thymine 
which exceeds the estimated electronic coupling between 
base pairs. In addition, DNA's environment, e.g. its dy- 
namical interaction with water and counterions, can have 
a substantial impact on its conductivity. From molecu- 
lar dynamics simulations, the average displacement of a 
base pair in DNA is about 0.3 — 0.4A [2] which is a tenth 
of the Watson-Crick spacing and an order of magnitude 
higher than in crystals at room temperature. All these 
properties make DNA a highly dynamic system and it 
remains unclear how well traditional concepts from solid 
sate physics can describe it. 

DNA has recently been the object of numerous single 
molecule studies to ascertain its conducting properties 



directly. Results have varied dramatically, from finding 
a fully localized insulator [3], a wide gap semiconductor 
with coherent extended states [4], a metal at room tem- 
perature with ohmic current voltage characteristics and a 
conductance comparable to doped polyacetylene[5], and, 
most surprisingly, a proximity effect superconductor at 
very low temperatures [6]. Recently, self-assembled net- 
works of DNA have been developed with quite reason- 
able conductance properties over distance scales of 50 
nm or less, and high conductance induced by chemical or 
gated doping [7, 8]; these networks can be used to make 
transistor devices. Indeed, a single molecule DNA FET 
(field-effect transistor) has recently been constructed [8]. 
Motivated by these wide ranging experimental results on 
the conductivity of DNA, we have embarked on a theo- 
retical effort to ascertain what conditions might induce 
such remarkable behavior. Our focus here is to exam- 
ine whether any likely DNA-structures or environments 
can yield reduced activation gaps to conduction or en- 
hanced electronic overlaps. In particular (sec. II), we 
have studied a hypothetical stretched ribbon structure, 
A, and B-form DNA, and the effects of counterions (espe- 
cially sodium and magnesium) with and without water. 

We have used a combination of fully ab initio density 
functional theory (DFT) code (SIESTA [9]), and a pa- 
rameterized Hiickel-Slater-Koster[10, 11] modeling using 
only atomic Pz orbitals [z the long axis of the molecule) 
of the bases. 

We have been unsuccessful in identifying structural 
motifs that might engender metallic conductance in 
DNA, but we give an explanation why DNA in the A- 
form provides less effective electronic coupling along the 
helical axis than B-DNA. While direct matrix elements 
between two highest occupied molecular orbitals (HO- 
MOs) or lowest unoccupied molecular orbitals (LUMOs) 
for model A-(B-)DNA are of order 0.01 eV (0.1 eV), indi- 
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vidual matrix elements between Pz orbitals on well con- 
tacted atoms in adjacent base pairs can generally exceed 
1 eV. 

We also demonstrates the possibility of small activa- 
tion gaps of order ksT with a first principles study of 4 
base pair long B-DNA. We discuss the cases of sodium 
and magnesium counterions, with (wet DNA) or without 
water (dry DNA). 

We have also studied the influence of longitudinal and 
torsional vibrations upon the electronic structure of DNA 
(sec. Ill), and estimated reorganization energies, Franck- 
Condon factors, and charge transfer rates between adja- 
cent bases (sec. IV). We find good agreement between 
our estimated rates and recent experimental data[12] as- 
suming that torsional (twist) vibrations limit the charge 
transfer most significantly. Finally, we conclude in sec- 
tion V. Some of the mathematical details can be found 
in appendix A. 



II. THEORETICAL SEARCH FOR HIGH 
CONDUCTANCE DNA 

A. Structural influence 

In this section we demonstrate that the electronic cou- 
pling between base pairs is highly sensitive to the DNA- 
structure and thus to the influence of its environment, 
which can modify the structure. 

In biological situations, there are several double helical 
conformations of DNA depending on the humidity and 
salt concentration. We consider the A- and B-forms. The 
B-form predominantly occurs in vivo, has a unit cell of 
10 base pairs (neglecting sequence effects), a helical rise 
of about 3.375A per base pair, and a twist angle between 
adjacent base pairs on average of 36.0° [13]. Since the 
base pairs are approximately perpendicular to the helical 
axis, the base pair separation is close to the helical rise. 
The A-form exists at low humidity in the presence of 
some salts. The double helix is relaxed to 11 base pairs in 
a unit cell exposing more of the hydrophobic core of base 
pairs and portions of the sugar units of the backbone. 
The helical rise and base pair separation are 2.56A and 
2.425A, respectively. The twist angle is only 32.7° [14]. 

There are also other known conformations. In order 
to obtain straight DNA molecules for conductivity mea- 
surements, one often uses a technique called molecular 
combing [15]. When a DNA molecule is anchored at one 
end on a hydrophobic surface like polystyrene, a receding 
meniscus can stretch DNA by a factor of 1.7, or if grafted 
at both ends to the surface even by a factor as large as 2.1. 
The capillary forces are typically > 160piV and can be- 
come as large as 500piV which results in breaking of the 
molecule. The stretching has been theoretically exam- 
ined with molecular dynamics simulations using classical 
force fields [16]. Since the phosphate-phosphate distance 
along the backbone is about 7 A and the base pair sep- 
aration in B-DNA is about 3.4A, stretching by a factor 




FIG. 1: Ideal 2-dimensional ribbon structure for DNA 
stretched by a factor of 2. The box indicates the unit cell 
of height 6. 75 A. There are only very few good electronic con- 
tacts. 



of 2 can, in principle, lead to a complete unwinding of 
the double helix. Since just stretching would weaken the 
inter base pair electronic coupling, we examine the possi- 
bility of base pairs which lie in the plane of the backbone 
strands (fig.l). In this case, the base pairs can come 
close again and the kinetic energy of the tt electrons can 
in principle be lowered if there was decent ppn overlap be- 
tween the Pz atomic orbitals at the closest contacts (with 
the z-axis now perpendicular to the ribbon). Ribbon-like 
structures are indeed found by K. M. Kosikov et al [16], 
e.g. the structures Ai or by using a classical force 
field. 



B. Ab initio estimates of electronic couplings 

In order to get a realistic picture of the electronic 
coupling strength between two base pairs, the small- 
est structural unit through which charge has to go, 
we use the DFT code SIESTA [9]. Since its aim is 
efficiency and large systems, it uses Troullier-Martins 
norm-conserving pseudo potentials [17] in the Kleinman- 
Bylander form[18], and a basis set of numerical atomic 
orbitals using the method by Sankey and Niklewski[19]. 
We used a double-^ basis set with polarization orbitals 
(DZP), the largest basis within SIESTA, and the gener- 
alized gradient approximation (GGA) for the exchange- 
correlation energy functional in the version by Perdew, 
Burke and Ernzerhof [20] . The structures for B and A- 
DNA were obtained from x-ray diffraction data ([13] and 
[14], respectively). Since the hydrogen atoms were miss- 
ing in the A-DNA structural data, they were added and 
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relaxed with the conjugate gradient method. The ideal 
2-diniensional ribbon structure of figure 1 was created 
artificially in order to get most effective ppir overlap and 
a small unit cell containing only a single base pair. As 
for the 2 base pair calculations, only methylated[21] base 
pairs without the backbone were used for simplicity. 

Consider two molecular orbitals (MOs), a from base 
pair 1, and /3 from base pair 2. (y , f3 are eigenstates from 
SIESTA for individual base pairs. They can be expanded 
in terms of the original atomic orbitals of SIESTA, e.g. 
\a >= X]iI°iC"|ao >i with Na being the dimension of 
the base pair 1, and cf and \ao >i the LCAO (linear 
combination of atomic orbitals) coefficients and atomic 
orbital basis states, respectively. 

The electronic coupling between states a and (3 can 
be computed as follows using formulas for non-orthogonal 
states [22], viz 

tafS = [i'„;3-0.5(C+i;3/3)5a/3](l-O"' W 

*W = E E '^^H^' S^0 = 12I1 (2) 

i=l j=l i=l j=l 

where 7 and S can take the values a or {3 according to eq. 
1. The total Hamiltonian H (overlap S) matrix of the two 
base pair system with dimensions (Ar„ + A^^) x (Na + Nj^) 
can be obtained from SIESTA. It can be thought of as 
consisting of two diagonal blocks iJ"", Hl^^^ (5"", S^'^) 
representing the intra-base pair couplings (overlaps) and 
one off-diagonal block H"^ {S"^) representing inter-base 
pair couplings (overlaps). 

The values for the electronic couplings taf3 between two 
HOMOs, as well as two LUMOs, are shown in table I. 
As for the B-form, they arc in general small and in good 
agreement with earlier ab initio studies [25-27], but for 
A and ribbon form they are much smaller than for the B- 
form. For the ribbon structure (including the backbone), 
we also carried out a band structure calculation after 
fully relaxing the geometry. The bands were calculated 
at 30 k points in the long direction using a Brillouin zone 
sampling at 3 k points. There was no sign of dispersion 
in agreement with the weak electronic coupling. That is, 
once stretched the base pair states localize. 

Thus, our conclusion is that dehydrating (A confor- 
mation) and stretching of DNA should lead to even more 
highly localized states and insulating behavior than in 
B-DNA. For the ribbon, the main reason is that there 
are only very few (« 3) close contacts (~ 2.7 A) between 
neighboring base pairs (sec also figure 1). For A-DNA, 
the weak coupling might seem a bit surprising since the 
base pair separation is much smaller than in B-DNA, but 
it can be understood in terms of less effective stacking. 
This is explored in the next subsection. 



C. Slater-Koster modeling of coupling strength 

In order to understand the results in the previous sub- 
section more quantitatively, we have performed model 



TABLE I: Transfer integrals in eV between HOMOs (tn) and 
LUMOs (ti) of two base pairs. The calculations are done with 
SIESTA and DZP basis set. The expression dimer always 
stands for a two base pair long DNA sequence, e.g. the GG 
dimer stands for 5'-GG-3' DNA. 
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tL 


tH 


tL 


GG 


0.0069 


-0.0006 


-0.1409 


0.0525 


GG(2.425A)(ref. [23]) 






-0.6922 


0.2548 


GG(0.0°)(ref. [24]) 






0.2385 


0.3233 


5'-AG-3' 


-0.0153 


0.0060 


-0.0710 


0.1124 


5'-GA-3' 


-0.0113 


-0.0010 


-0.1871 


0.0472 


AA 


0.0310 


-0.0120 


-0.0695 


0.1054 


Ribbon 


GG 


0.0039 


0.0083 







calculations varying the angular amd spatial separation 
of bases. These are more readily performed than full 
ab initio calculations and often provide more physical 
insight. We restrict ourself to just atomic Pz orbitals, 
where the z axis can be defined locally by the normal 
of the corresponding base. In this case, two Pz orbitals 
from different base pairs couple by ppcr and pp7r interac- 
tions, which we model with semi-empirical Slater-Koster 
theory[10, 11] 



ppX — VppX 



rrii 



(3) 



We note rjppa- > and ?7pp7r < 0. Hero, d and m are dis- 
tance and electron mass, and {h'^)/{mdl) = 7.62eV. The 
exponential distance cut-off Rc and r] are parameters to 
be determined by matching to ab initio results as dis- 
cussed below. The interatomic matrix element between 
two "parallel" Pz orbitals depicted in figure 2 is 



sin (j>Vppa 
md'i 



f cos (t>VppTr 

ijlppa + |??pp7r|) 



(4) 



P + z"^ 



I ^ppTT I 



According to this formula, the combination of close base 
pair separation and poor contacts reduces the electronic 
coupling between base pairs, i.e. can lead to cancellation 

effects. 

A general formula for non-parallel orbitals, which is 
used for all our following model calculations can be found 

in the Appendix, cq. (A3). The normal of each base 
is taken for the direction for all pz orbitals within that 
base, and is obtained by fitting a plane to each single 
base whose coordinates were obtained from x-ray diffrac- 
tion data for A and B-DNA. This formula constitutes a 
straightforward way to calculate inter - base pair cou- 
plings for a variety of spatial and angular separations. 

The Slater-Koster parameters rjpp^, ripper and cut-off 
Rc were determined from fitting to DFT data as follows: 
Starting from the DZP basis of SIESTA, we first iden- 
tiiy the 2pz orbitals. These are the first-C orbitals of 
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FIG. 2: Iritor-atomic matrix elements from Slater-Koster the- 
ory: two atomic pz orbitals on adjacent base pairs couple by 
positive ppa and negative pp7r interactions. 



Pz symmetry. (The sccond-C orbitals correspond to ex- 
cited atomic states and have rather small contributions 
to the wavefunctions and matrix elements [< 20%].) We 
consider the reduced SIESTA Hamiltonian matrix in this 
basis, wliicli has two diagonal blocks representing the on- 
site and intra-base pair couplings and one off-diagonal 
block representing inter-base pair couplings. The Slater- 
Koster matrix elements calculated with eq. (A3) are di- 
rectly fitted to the off-diagonal block matrix elements. 
The fitting was done with simulated annealing, which 
gave better (and still physical) results than e.g. with the 
Powell algorithm[28]. Finally, since the DFT code uses a 
global z-axis along the helix which is different from the lo- 
cal z-axis due to an inclination and a propeller twist angle 
associated with each base, we turn to parallel reference 
base pairs for the fitting process. In the case of A-DNA 
we use two parallel G ■ C base pairs with a separation of 
2A3A, while for B-DNA ones with a separation of 3.34A. 
The fitted parameters are ripp„ = 2.93, ripper — —0.73, 
Rc = 1.161 for A-DNA, and rjppa = 5.27, rjpp^ = -2.26, 
Rc = 0.871 for B-DNA. 

In order to further stress the difference between A and 
B-DNA, figure 3 shows a distribution of all possible Pz 
interatomic matrix elements between two G ■ G base pairs 
in the A-form and in the B-form. Since guanine and 
cytosine have 11 and 8 Pz orbitals, respectively there are 
a total of 19 X 19 = 361 matrix elements. (Adenine and 
thymine have 10 and 8 Pz orbitals respectively.) Figure 
3 shows that there are more good contacts, e.g. twice 
as many above 0.75 eV, in the B-form than in A-form. 
Furthermore, there is a shift to (negative) ppw interaction 



40 

30 

20 - 



g 10 

o 

I 



40 



30 



20 



10 



Slater-Koster matrix elements Vy[eV] 



ADNA 



a-type 




^ - rH r*""^ n n ^.n i-j n 



0.25 



0.5 



0.75 



1.25 



BDNA 



1 a-type 



0.25 



0.5 



0.75 1 



1.25 



FIG. 3: Distribution of Slater-Koster inter-atomic matrix el- 
ements between pz orbitals of two stacked G ■ C base pairs. 

There arc a total of 361 matrix elements. The smallest peak 
height corresponds to a single matrix element. Details about 
the distributions, mean: A 0.09, B 0.13; standard deviation: 
A 0.21, B 0.23; number of matrbc elements above 0.75eV: A 
6, B 13. 



in the A-form, although the single largest matrix element 
(« 1.3 eV) arises in the A-form. This is presumably due 
to a single optimal contact at the shorter distance. 

Table I shows the electronic couplings between specific 
MOs which correspond to HOMOs and LUMOs. A way 
to understand the smallness of Ihomo {in) and tLUMO 
{tL) for A-DNA relative to B-DNA is to note that the 
distribution (fig. 3) is clearly more centered around zero 
in case of A-DNA, while for B-DNA it is more asymmet- 
ric with mostly positive matrix elements. The electronic 
coupling is a weighted average of these matrix elements 
with the weights being the product of two LCAO coef- 
ficients. It can be expected that the difference in the 
distributions leads to small electronic couplings for A- 
DNA. 

Since A-DNA appears clearly insulating, as already in- 
dicated by the small bandwidth reported in ref. [3], we 
restrict ourselves to B-DNA in the subsequent sections. 



D. Effect of counter ions ions and water on 
electronic structure 

DNA is a highly charged molecule which can only be 
stable if the negatively charged phosphate groups are 
neutralized by positive counterions or by polarized water 
molecules from a buffer solution. In order to measure its 
conductivity DNA has initially to be dried, normally by a 
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flowing N2 gas [7]. This can result in as few as 2 to 3 wa- 
ter molecules per nucleotide[29], which is not even high 
enough humidity to obtain the A-structure. In principle 
dried DNA can turn into a disordered structure [30] prob- 
ably not suitable for charge transport through the base 
pairs at all. Evidence for that is found by STM images of 
dried DNA on mica, which shows a base pair separation 
of ^ 7.2A [31], and also by AFM images on Si02 which 
show a DNA sample height of only 0.5 nm instead of the 
diameter 2nm of the usual double helix[32]. 

Different research groups in general work with dif- 
ferent buffer solutions: activated (0.12 - 0.2 cV) hop- 
ping conductivity with an initial sodium ion buffer was 
observed by Kawai's group [8] using /im long poly(G)- 
poly(C) /poly(A)-poly(T) DNA, and by Rakitin et al [33], 
as well as Tran et al [29] using A-DNA. On the other 
hand, proximity effect superconductivity was observed 
by Kasumov et al [6], who used a few 16-/.tm-long A-DNA 
molecules taken from a magnesium buffer solution. 

Motivated by this, we have examined the effects of var- 
ious counterions (protons, Na, Mg) and water on the 
electronic orbitals of DNA. We use a 4 base pair long 
5'-GAAT-3' structure obtained from classical molecular 
dynamics (MD) with the AMBER 4.1 force field. The 
structure was taken from a snap shot of the Dickerson 
dodecamer with thousands of water molecules [2] and 
was provided to us by S. Dixit. For wet DNA we keep 
the first and second solvation shell (~ 25 water molecules 
per nucleotide) , which corresponds to keeping all the wa- 
ter molecules within a 4.5 A radius of the DNA atoms. 
Although B (wet) DNA exists only above 13 to 18 wa- 
ter molecules per nucleotide, we also examine completely 
dry DNA by removing all the water molecules but keep- 
ing the same B-DNA-structure. From our own experi- 
ence with DFT calculations, DNA orbitals are generally 
rather localized on various groups like phosphate, sugar, 
counterions, and bases. Thus it might still be possible 
that this artificial dry DNA shows the correct order of 
energy levels and hence gives valuable information about 
activation gaps. The location of the counterions are kept 
the same as for wet DNA, i.e. some are further outside 
the DNA helix while others are at the phosphates or in 
the grooves. Counterion distributions are discussed in 
ref. [34]. For comparison we also considered protonated 
DNA of the same structure, since this is what theorists 
usually use for neutralized DNA. This has the big ad- 
vantage of keeping the number of atoms small. The pro- 
ton counterions are placed at the phosphates where they 
bind covalently. All other metal counterions and water 
molecules are removed. 

We optimized the geometries with a few conjugate gra- 
dient steps (time consuming) , but the forces still large (of 
order eV/A). On can argue that this out-of-equilibrium 
structure might nevertheless resemble DNA in a solution 
at finite temperature where forces are always non-zero 
due to the motion of atoms and molecules. As a further 
check for dry DNA, we placed the metal ions at different 
locations and observed no principle changes occurred in 
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FIG. 4: PDOS of protonated 4 base pair long B-DNA; black: 
total DOS, red: pz orbitals of atoms in base pairs, green: 
phosphate, blue: sugar. The Fermi energy is set to eV. 

the band structure or the MuUiken population [35] of the 
ions. For DNA we use a double-C basis set except for 
phosphorus, hydrogen involved in hydrogen bonding and 
the counterions for which we also include polarization or- 
bitals. For the water molecules (more than 200) we only 
use a minimal, single-C (SZ) basis. The density of states 
(PDOS) is then projected onto atomic orbitals for ener- 
gies around the Fermi energy ep [36] to see which atomic 
orbitals can contribute to charge transport and what ef- 
fects the counterions have on the eigenvalue spectrum. 

The PDOS of the ith atomic basis orbital projected on 
the nth MO is 

Nb 

p7i^) = Y.c: S,, 6{e-e„), (5) 
j 

where Nb is the dimension of the basis set, cf,c" are 
the LCAO coefficients and Sij are the atomic overlap 
matrix elements. We project on the atomic pz orbitals 
of the elements C, N and O of all the bases (resulting in 
TT and TT* MOs), all orbitals of the five phosphate atoms, 
of the sugar atoms, counterions, and also whole water 
molecules. 

The result for the protonated DNA is shown in figure 
4. The atomic Pz orbitals (red) of the bases form an 
expected tt - tt* gap of about 1.9 eV around ep (set to 
zero) . One has to keep in mind that DFT usually under- 
estimates the gap in insulators by about 50 to 100%. The 
HOMO, mainly guanine, is about 0.4 eV higher than the 
next lower occupied MOs which are mostly adenine. The 
LUMO (cytosine) is about 0.2 eV lower than the next 
higher unoccupied MO (thymine). The occupied (unoc- 
cupied) phosphate states (green) and sugar states (blue) 
on the other hand are energetically lower (higher) than 
the TT (tt*) orbitals. 

The results for sodium are shown in figure 5. Parts a) 
and b) contain the PDOS of wet DNA, i.e. part a) shows 
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FIG. 5: Effects of sodium counterions and water on the molec- 
ular orbitals of 4 base pair long B-DNA. Part a) and b) show 
the PDOS of wet DNA; a) red: Pz, green: phosphate, blue: 
sugar, yellow: sodium; b) blue: all water molecules, green: 
water molecules near (< 3.5^4) phosphates, yellow: water 
molecules near sodium, red: water molecules near bases. In 
part c) all the water are removed (dry DNA). 



the projection onto the bases (pz), phosphates, sugars, 
and sodium orbitals, while b) shows the projection onto 
water. In a) one can see a tt - tt* gap of about 2.2 eV. It 
has slightly increased in the presence of water compared 
to protonated DNA. The highest tt states around -1.6 eV 
(relative to e^) are mainly guanine, the next lower group 
of occupied states —1.9 eV) are quite delocalized on 
several adenines with some guanine admixture. The first 
TT* state has mainly thymine character, the next higher 
also a little bit of cytosine. In the presence of water the 
TT and TT* MOs seem to be more extended than in case of 
protonated DNA, where the MOs are more localized on 
individual bases. 

Although the tt - tt* gap is rather large, the real HOMO 
- LUMO gap is surprisingly small (54 meV). Electrons 
could possibly be excited from water levels below e p into 
unoccupied water and sodium states. The presence of wa- 
ter and counterions might result in a weak conductance 
if a voltage is applied. Part b) of figure 5 separates the 
water PDOS (blue) into contributions coming from wa- 
ter molecules near the phosphates (green), near sodium 
(yellow), and near the bases (red). As one can see, the 
water states around ep are mostly from water molecules 
near the sodium ions. The LUMO is a mixture of water 
and sodium. Furthermore, the negatively charged phos- 
phates and positively charged sodiums are affecting the 
electronic structure of water drastically, since the several- 
eV large insulating gap of water is reduced to about 1 



eV (0.7 eV in case of magnesium). A MuUiken popu- 
lation analysis of sodium reveals that sodium is on av- 
erage a -1-0.88 charged ion. This is consistent with first 
principle studies of solvated metal ions, which indicate 
that solvated sodium is nearly a -fl-ion almost indepen- 
dent of the water coordination number [3 7]. Regarding 
the sodium distribution, two sodiums are several A out- 
side the DNA double helix, but the MuUiken population 
is hardly any different for these. 

In order to see the effects of water on DNA more dras- 
tically, part c) of figure 5 shows the DOS of DNA with 
all the water molecules removed. Since the sodium ions 
(now only -fO.66) are not directly located at the phos- 
phates, they cannot screen their negative charges. Re- 
pulsive interactions likely between the oxygens of each 
phosphate group increase their energy (green) even be- 
yond the TT states (red). Although the tt - tt* gap is 
about 2 eV, there are sodium states only 39 meV above 
the HOMO (phosphate). The PDOS is remarkably re- 
producible for different counterion configurations. One 
can speculate that due to the small electron excitation 
gap hole hopping through the backbone might be possi- 
ble. Even if the B-DNA-structure is unstable under per- 
fectly dry conditions, the quasi one-dimensional pathway 
for conduction through the backbone might still survive. 
One the other hand, the insulating sugars (blue) between 
the phosphates do not contribute much to the relevant 
states and constitute tunneling barriers. 

For magnesium counterions, the result is shown in fig- 
ure 6. Parts a) and b) contain the PDOS of wet DNA. 
Part a) shows the projection on DNA, while b) shows the 
projection onto water. In a) one can see a tt - tt* gap of 
about 2.2 eV which is essentially the same as in the case 
of sodium. The highest tt states at about -2.2 eV are a 
mixture of guanine and adenines, and 0.8 eV above an 
occupied adenine/ thymine MO. The tt - it* gap is again 
much bigger than the real HOMO - LUMO gap, which is 
only 62 meV. Magnesium and water states below ep could 
function as electron donors with electrons being excited 
into unoccupied tt* MOs leading to a possible electron 
conduction mechanism. From part b) of figure 6 one re- 
alizes however, that the water states right below ep are 
not from water molecules in close proximity to the bases 
making the electron donation harder. The occupied wa- 
ter states from water near the bases are about 1 eV below 
the TT* MOs. A MuUiken population analysis of magne- 
sium reveals that magnesium is not a -|-2-ion, but has a 
reduced -1-1.25 charge. One magnesium is more distant 
from the DNA molecule and is completely surrounded by 
water (7 water molecules within a radius of 3A). It has 
even a smaller (+1.0) charge. The smaller charge (less 
than -t-2) is again consistent with first principle studies 
of solvated metal ions, which indicate that magnesium 
solvated by 4 water molecules has a -1-1.1 charge, but less 
solvated magnesium has a larger positive charge closer to 
+1.7 [37]. 

Part c) of figure 6 shows the DOS of DNA with all 
the water molecules removed. Although the tt - tt* gap 
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energy [eV] 

FIG. 6: Same as figure 5 except for magnesium counterions 
instead of sodium. Part c) : there is a single magnesium state 
right above the Fermi energy {ep = 0). 



is about 2.8 eV, there is a single magnesium state (not 
from the further outside one) only 11 meV above the 
HOMO (phosphate). The PDOS is again remarkably re- 
producible for different counterion configurations. 

Finally, we try to give a possible explanation for the 
accumulation or pinning of phosphate states right be- 
low the metal states for dry DNA. Left unscreened, the 
occupied phosphate states will rise in energy due to re- 
pulsive interactions between negatively charged oxygens 
of a phosphate group. However, they cannot get higher 
than all the empty metal states. This is because if they 
were higher, the phosphates would dump electrons in the 
metal states and reduce their repulsive interactions. This 
would lower their energy, but if below the metal states, 
the phosphates would accept the electrons again, and the 
same process would repeat itself leading to the pinning- 
effect. 



III. DYNAMICAL INFLUENCE ON 
ELECTRONIC STRUCTURE 

At finite temperature the base pairs of DNA will oscil- 
late about their equilibrium positions. Within our model 
from section II C we are ultimately interested in the tem- 
perature dependence of the n - n* band gap 



A 



DNA 



(T) 



A 



single bp 



2(<iHOMO> (T)+ <twMO> (T)), (6) 



with Asingio bp being the temperature independent 
HOMO-LUMO gap of an isolated base pair, and <t> (T) 



describing the temperature dependent average electronic 
coupling between adjacent base pairs. (In this model 
calculation, the tt - tt* gap corresponds to the HOMO- 
LUMO gap, because we do not consider water here.) 
The goal is to understand the experiment by Porath et 
al who measured a strong increase of the voltage gap 
Uc (the minimal bias voltage above which conduction 
sets in) with temperature (2-4 eV per 300 K)[4]. In 
the experiment, short (30 base pair long) homogeneous 
poly(G)-poly(C) DNA molecules showed band-like semi- 
conducting behavior. In our calculation, we assume that 
displacements along the helical axis and the twisting mo- 
tion are the most important degrees of freedom which 
affect the electronic coupling. For small changes in the 
variables, one can assume that these two degrees of free- 
dom are independent. 

Let (f> and u describe the deviations from equilibrium 
twist angle and base pair separation. Within our Koster- 
Slater-model in the CNDO (complete neglect of differ- 
ential overlap) approximation, the transfer integral be- 
tween two MOs of adjacent base pairs is given by eq. 1 
with 5i2 = 0. (From our first principles calculations, 
0.001 < Si2 < 0.01.) Thus we have ti2 = Y.^,j Hlfc\c] 
where the LCAO coefficients are summed over each of the 
two base pairs and kept fixed for all cj) and u. They are 
given by the first-C orbitals of symmetry from a single 
base pair DZP calculation, and are corrected by a factor 
cos""'^ 9 where 9 — 8.0° is the angle between the local 
base normal and the helical axis. As noted earlier, the 
second-C orbitals yield < 20% corrections to this. The 
matrix H^^ is obtained as described in sec. II C and de- 
pends on (j) and u through the Slater-Koster form of the 
matrix elements. The tilde indicates that the reduced 
(first-C) basis set was used. Since we calculate transfer 
integrals in CNDO, we do not need on-site energy ma- 
trices H^^^ , H^^ here. Nevertheless, these can be easily 
obtained in the same fashion. The resulting model can be 
used for very large scale quasi-self-consistent-field quan- 
tum calculations of long DNA sequences - similar to the 
one used for electron transfer rate calculations between 
various donors and acceptors in DNA[38]. 

In figures 7 and 8 we show the electronic coupling be- 
tween the frontier orbitals as a function the variables 4> 
and u. One of them is held fixed at the equilibrium value, 
while the other one is varied. Consider first the variation 
with the angle </> (fig. 7). At </> = —36°, i.e. when the 
two base pairs are aligned perfectly parallel, the coupling 
is maximal and positive. Let us focus on a GG dimer. 
At = —36°, there are 19 optimal contacts right on 
top of each other, and the transfer integral is approxi- 
mately given hy t ^ Yli {ciY {ci := cj = cf) where 
Hf^ > is solely due to ppa interaction. Note the inter- 
esting property that there are sign changes of t((/i) with 
(f>. These provide a direct explanation for the local min- 
ima of \t\ observed in earlier ab initio studies [25, 26]. In 
B-DNA, essentially all the interaction matrix elements 
{Hl^ > 0) are positive {ppa dominated), as can be seen 
in figure 3 and this is true for arbitrary twist angles (j). 
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twist angle [deg] 

FIG. 7: Change of electronic coupling between frontier or- 

bitals of two base pairs with relative twist angle 4> about the 
equilibrium position. The base pair separation is kept fixed 
at = 3.3751. 
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FIG. 8: Change of electronic coupling with relative displace- 
ment u about the equilibrium position. The twist angle is 
kept fixed at 36.0°. 



Since for larger twist angles, i.e. ~ 0, the transfer in- 
tegral is approximately given by i ~ X^i/j ^Ij^l^j^ 
sign change of t can only come from the signs of the 
LCAO coefficients (phases of the Pz orbitals). This can 
lead to positive or negative t, or to a complete cancella- 
tion. 

Coming back to the temperature dependence of t and 
the band gap, we assume classical harmonic potentials 
with spring constant K and shear constant S for each 
couple of base pairs. They can be estimated using the 
equi-partition theorem from classical statistical mechan- 
ics and estimates for the standard deviation from molec- 



< > -- 



ular dynamics [2], 

0.375^1^ 

— > K = 0.184 

— > S 



K 



eV 
l2 



<02 >= 



294.8 
knT 



pN 



4.6-10- 



S 

eV 



deg 



,2 ■ 



(7) 

(8) 

(9) 
(10) 



With these estimates for K and S the average electronic 
coupling can be calculated 



<t>{T) 



(11) 



with /3 = I/UbT. The transfer integral only depends 
on the differences u = U2 — u\ and (f) = (f>2 — <pi, where 
Ui and arc the displacement and twist angle of the 
ith base pair. The results are shown in figures 9 and 
10. Two things are remarkable about this result. First, 
the temperature dependence is very weak (a few meV 
change over hundreds of K). Second, harmonic displace- 
ment along the helical axis alone can lead to an increase 
in magnitude of the electronic coupling with increasing 
temperature T, while twisting motion on the other hand 
has the opposite effect leading to a cancellation. 

This can be understood by looking at figures 7 and 
8. Since t{(j) = 0) is closer to a local maximum (at least 
for HOMO states), a broadening of the gaussian distri- 
bution e~^^'^ 1"^ with increasing temperature about 
(/) = will give more contribution from smaller t. On 
the other hand, a broadening of the distribution for u in 
the form e~^^" 1"^ will lead to more contributions from t 
with larger magnitude, since t{u) increases exponentially 
in magnitude when lowering u. 

The weak temperature dependence of the electronic 
coupling, and hence of the band gap from eq. (6), are 
in strong contrast to the experiment by Porath and co- 
workers. It also calls to question the earlier analysis by 
Y. Berlin et al [39] who tried to explain the temperature 
dependence of the observed voltage gap Uc from a reduc- 
tion of the bandwidth by only including torsional motion 
of base pairs. Our results imply a temperature indepen- 
dent Uc due to the cancellation of the two modes. The 
extremely strong increase of the gap with temperature 
could well be due to other effects [40]. The smallness of t 
does not allow much decrease of the bandwidth anyway. 



IV. ESTIMATION OF PARAMETERS FOR 
CHARGE TRANSFER 

Single electron or hole transfer is likely to be very sen- 
sitive to the actual motion of base pairs. Experimental 
evidence comes from the observation of two different time 
scales (5 and 75 ps) in charge transfer experiments[41]. 



9 



0.04 - 

LUMO GC - LUMO GC 

0.035 - 



— K+S 




0.025 - 



I I \ I \ I \ I \ I \ I I 

50 100 150 200 250 300 

T[K] 

FIG. 9: Boltzmann averaged electronic coupling between LU- 
MOs of two G ■ C base pairs as a function of temperature T. 
The effects of a spring constant K and a shear constant S are 
shown. 



TABLE II: Inner reorganization energies Ex^i, moments of in- 
ertia / and nuclear frequencies Vn of single base pairs. The 
calculations are done with SIESTA and DZP basis set. E\^i is 
determined by the energy difference of a hole/electron doped 
base pair between doped and undoped relaxed geometries. 
The moment of inertia is estimated with / = X] "^i'"? where 
the sum is over all atoms of a base pair and are mea- 
sured from the helical axis of B-DNA. The nuclear frequencies 
are derived with E\^i = 2{ni')^ ^ A*i,i+iWi,i+i with and 
w<,i+i being the reduced mass and change in bond length of 
atoms i and neighbor i+l when comparing doped with un- 
doped relaxed structures. 





GC 


A-T 


Ex,i [eV] 


hole 


0.18 


0.08 




electron 


0.28 


0.20 


J[10-^^kgm^] 


5.0 


6.2 


i/„[10"l/s] 


hole 


3.97 


3.64 




electron 


5.09 


2.74 
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> 
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-0.15 
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Konly 
S only 



■ 50 100 150 200 250 300 

T[K] 

FIG. 10: Boltzmann averaged electronic coupling between 
HOMOs of two G • G base pairs as a function of temperature 
T. The effects of a spring constant K and a shear constant S 
are shown. 



The longer time stems presumably from a necessary re- 
orientation of base pairs in order to make charge transfer 
possible[41, 42]. Theoretically, torsional acoustic modes 
are indeed found to be extremely soft (< 20cm^^) [43] 
and hence slow, which is consistent with dynamic Stoke 
shifts in fluorescence spectra[44]. According to figure 7, 
the magnitude of the transfer integral oscillates between 
approximately and 0.1-0.2 eV when one assumes a stan- 
dard deviation of the twist angle of about 10° (actually 
more like 7.5"). Under this assumption electron tunnel- 
ing is possible on a time scale of the oscillation time T5 
which can be estimated as Tg = I-k^I jS^ where S is 
the shear constant and / = the reduced 

moment of inertia of two base pairs. 

More specifically, single charge transfer rates between 
two base pairs can be estimated from the classical Marcus 



formula [45] (neglecting nuclear tunneling effects) 

kcT = yn^de~'^^^\ (12) 

where and b.E* = {Ex + Eof/{'iEx) are the effec- 
tive nuclear frequency and energy to reach the transition 
state, while k^i describes the electronic transmission coef- 
ficient. E\ is the total reorganization energy and Eo <Q 
is the reaction energy which is zero for charge transfer 
between homogeneous base pairs. 

Wc expect the charge transfer between neighboring 
base pairs to be adiabatic, i.e. Kei ~ 1, since the Landau- 
Zener probability P for crossing the surfaces at the tran- 
sition state is close to unity with [46] 



and 



P 

1 



exp{—2n'y) 



ShunV-^kBTAE* ■ 



(13) 

(14) 



This is because 7 is large due to the relatively large trans- 
fer integral t. Inner shell reorganization energies and 
nuclear frequencies associated with adjusting the bonds 
between atoms of a single base pair can be calculated by 
the method of A. Klimkans and S. Larsson described in 
[47]. The results are shown in table II, which also con- 
tains estimates for moments of inertia. The frequencies 
are quite similar to in-plane E2g frequencies (4.8 lO^^Hz) 
in graphite [48], and our calculated inner shell reorgani- 
zation energies lie between those of benzene (0.26 eV) 
and anthracene (0.07 eV) [47], as one would expect from 
the number of aromatic rings involved. Since the total 
reorganization energy is estimated from experiments to 
be only about 0.4 eV [49] which is close to our values for 
the inner reorganization energies from table III, it sug- 
gests that the base pairs must be well protected by the 
double helix formation and solvation effects are minor. 
The small reorganization energies also contribute to the 
adiabaticity by increasing 7 in equation 14. 
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TABLE III: Oscillation time of twisting motion Tg, total re- 
organization energy Ex ~Y1 Franck-Condon factor FC, 
and upper/lower limits of adiabatic charge transfer rates be- 
tween two neighboring base pairs (bps). The rates are cal- 
culated with eq.(12) and data from table (II) in the normal 
regime. The upper limit corresponds to a nuclear frequency 
1/v — l/vn + l/v^, the lower limit to v = 1/Ts- 



sequence(ref. [50]) 


GG 


AG(ref. [51]) 


AA 


Ts[ps] 


2.02 


2.13 


2.25 


Ex[eY\ 


hole 


0.36 


0.26 


0.16 




electron 


0.57 


0.48 


0.39 


FC (at 300K) 


hole 


0.03 


0.88 


0.22 




electron 


0.004 


0.01 


0.02 


limit 




khou [1/s] 


upper 


6.010'' 


1.710'-' 


4.1 10'^ 




lower 


1.5 10'" 


4.1 10" 


9.9 10'" 


^electron [1/^] 


upper 


1.1 10" 


1.810" 


3.1 10" 




lower 


2.1 10" 


4.7 10" 


1.010'" 



Upper and lower bounds for the charge transfer rates 
between two base pairs are shown in table 111. For the 
upper bound, we take f in equation (12) to be the inner 
base pair nuclear frequency, while for a lower bound we 
take it to be the base pair twist frequency. This slow pro- 
cess can be viewed as reaching an encounter complex with 
a frequency z/ = 1/Ts- We note the large Franck-Condon 
factor in table III for hole transfer between A ■ T and 
G-C (dimer AG). This is because this case is close to the 
so-called optimal electron transfer regime in which the 
reorganization energy (0.26 eV) is largely compensated 
by the reaction energy (0.20 ± 0.05 eV) [52] as estimated 
from analyzing electron transfer yield data. Note that 
this reaction energy is smaller than the difference of the 
bare oxidation potentials (0.4 — 0.5 cV) between isolated 
adenine and guanine bases in a polar solvent [53] , in gas 
phase [54] , or from ab initio calculations and Koopman's 
theorem[25, 55]. Furthermore, hole transfer is in most 
cases at least one order of magnitude faster and hence 
more efficient than electron transfer. 

In orck^r to compare with experiment turn to (hole) 
transfer between initially excited 2-aminopurine (an iso- 
mer of adenine) and guanine [12]. The measured charge 
injection time of lOps (corresponding to a rate of 10"'^"'^ 1 /s) 
is surprisingly close to our lower bound value 4.1 10^^1/s. 
Hence we speculate that the twisting motion of base pairs 
is rate limiting. 

Another simple application of the 2 base pair charge 
transfer rates of table III is the hole hopping between ho- 
mogeneous sequences, e.g. (A ■ T)jv, when one assumes 
(biased) random walk as the underlying mechanism of the 
charge migration. In this case, the rate shows a weak al- 
gebraic distance dependence [56] kcT ~ kiiop^~^, where 
N is the number intervening base pairs and 1 < r] < 2. 
Absolute rates can be predicted with our (lower bound) 
value fcnop = 9.9 10^°l/s w 10^^1/s. In order to compare 
with experiment, e.g. [57], one would have to include the 



initial transfer from the donor to the bridge, as well as 
the final transfer from the bridge to the acceptor which 
leads to a further reduction of the rate. General formulas 
are provided in refs. [52, 56]. 

Finally, since charge transfer between adjacent base 
pairs is likely to be adiabatic, one should avoid includ- 
ing these when fitting to a diabatic charge transfer rate 
kcT ^ exp(— /?i?). Including the adjacent base pair 
charge transfer in the fit leads to an under-estimation 
of /3. 



V. CONCLUSIONS 

In summary, we have studied the electronic proper- 
ties of DNA from several perspectives. Our initial mo- 
tivation was to come up with possible conditions which 
enhance the conductance of DNA. For example, we ex- 
pected a stretched ribbon-like DNA to be a potential 
candidate for a molecular wire. Unfortunately, that has 
not turned out to be the case. Indeed, in the ribbon-like 
DNA, the base pair separation can become quite small. 
However, the pp7r interaction appears not very effective 
as there are only very few good contacts. Using the ab 
initio code SIESTA we compared transfer integrals be- 
tween base pairs of A, B and stretched DNA. Our cal- 
culations suggest that A-DNA and stretched, ribbon-like 
DNA should support electrical current even less than B- 
DNA. 

On the other hand, a first principle study of 4 base 
pair long DNA in the B-form with different counterions 
(sodium, magnesium) under dry and wet 200 water 
molecules) conditions shows the possibility of small ac- 
tivation gaps of order ksT. Although the tt - tt* gap is 
several eV large, a small gap is formed mainly by water 
and counterion states in the case of wet DNA, or by back- 
bone and counterion states in the case of dry DNA. In 
the latter case, the backbone states rise in energy and 
lie higher than the occupied tt MOs due to repulsive 
interactions of the negatively charged backbone. Since 
small activation gaps (^ O.leV) are found experimen- 
tally for bundles [7, 33] or supercoiled DNA [8, 29] from 
a sodium buffer, this could be caused by residual water 
and condensed counterions providing an alternative path- 
way for charge hopping. To our knowledge, there is less 
experimental conductivity data available for DNA from 
a magnesium buffer. The most prominent case is the 
observed induced superconductivity [6] requiring truly 
extended metallic-like states. From our first principle 
studies of DNA with ions and water, we do notice some 
differences between magnesium and sodium ions. For 
wet DNA, electron excitation from occupied water and 
magnesium states into empty tt* states seems possible 
making electron conduction feasible. In addition, for dry 
DNA with magnesium ions, the activation gap between 
occupied backbone and unoccupied magnesium states is 
vanishingly small. In addition, there is a large DOS of 
phosphate states right below the Fermi energy available. 
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This pinning-efFect opens the possibility of a pathway for 
hole hopping through the backbone. 

In order to get a better microscopic understanding of 
the electronic couphng between adjacent base pairs and 
how it is affected by twisting and displacement motion, 
wc developed a Hiickel-Slater-Koster model with param- 
eters obtained from ab initio calculations. We find that 
the transfer integral is constrained by a competition be- 
tween ppcr and ppTT interaction and interference by phase 
factors of atomic Pz orbitals. The latter effect also leads 
to sign changes in the transfer integrals as a function of 
the torsional angle not anticipated in previous studies. 

Furthermore, the temperature dependence of the av- 
erage transfer integral and band gap is expected to be 
weak: First, the electronic coupling is already weak any- 
way, and cannot be reduced much. Second, twisting and 
displacement motion of base pairs have cancelling contri- 
butions. 

Next, we showed that charge transfer rates between ad- 
jacent base pairs are likely in the adiabatic limit. Thus, 
if one uses (modified) bases as donor and acceptor in 
charge transfer experiments, it is important to take this 
into account. When obtaining the decay length /3 of di- 
abatic charge transfer from a fit to an experimental rate 
as a function of distance, the data point corresponding 
to the donor and acceptor next to each other should not 
be included in the fit. 

Wc also calculated Franck-Condon factors, reorga- 
nization energies and nuclear frequencies to obtain 
absolute rate estimates. Furthermore, we suggest that 
the twisting motion might be rate limiting for hole 
transfer. Due to the adiabaticity of the two base pair 
rates, DNA might be a good medium for diffusive charge 
transport. 
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APPENDIX A: APPENDIX 

In the following we derive a general formula for the 
calculation of interaction matrix elements between two 
atomic pz orbitals belonging to adjacent base pairs. In 
general, the orbitals point in different directions, since 
the local z-axis is different for each base mainly due to a 
propeller twist angle, in case of A-DNA also because of a 
large inclination angle. We first expand the orbitals p^ , p'^ 
in a global coordinate system with the z-axis pointing in 
the direction of the helix 

Pz = apx + (ipy + "fPz (Al) 

P'z = Ol'Px + P'Py + l'Pz, 

where the expansion coefficients are the direction cosines. 
If (/, m, n) is the unit vector pointing from one orbital to 
the other, the Slater-Koster relations [10] 

= l''Vpp„ + {1 - l'')Vpp^ (A2) 
Eyy = m^Vppa + (1 - m^WppT, 
Ezz = n^Vppa + (1 - r?)Vpp-K 

Exz — ^'^{y^ppa ^ppn) 
Eyz — VriTliVppfj ^^p7r) 

can be used to obtain the following formula for the inter- 
action matrix element 

V = aa'E^x+PP'Eyy+jj'Ezz 

+{a/3' + /3a')Exy + {aj' + ja')Exz 

+ lP')Eyz. (A3) 

The Slater-Koster matrix elements Vppa and V^j,^ are 
given by eq. (3) and contain the distance dependence. 
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